Universality of striped morphologies 
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We present a method for predicting the low-temperature behavior of spherical and Ising spin mod- 
els with isotropic potentials. For the spherical model the characteristic length scales of the ground 
states are exactly determined but the morphology is shown to be degenerate with checkerboard 
patterns, stripes and more complex morphologies having identical energy. For the Ising models we 
show that the discretization breaks the degeneracy causing striped morphologies to be energetically 
favored and therefore they arise universally as ground states to potentials whose Hankel transforms 
have nontrivial minima. 



The study of pattern formation in simple systems has 
received much attention the last 20 years. Not only are 
the pictures visually arresting, producing reviews with 
considerable artistic qualities [l|, [2[, but spatially inho- 
mogeneous phases present difficulties for standard the- 
oretical methods, calling for new principles to describe 
the physics of systems exhibiting them |3[. Whatever 
weight can be assigned to the former as an explanation 
for the interest, the latter is justification enough, espe- 
cially as the experimental evidence of stripes, spots and 
checkerboards in practically important materials abound. 
Examples include lipid monolayers [4|], adsorbates on 
metals [Hj], and various magnetic fluids p, 0|. Striped 
phases are also hypothesized to play a role in the high- 
temperature superconductivity of transition metal ox- 
ides (HQ. 

Many experimental systems displaying heterogeneous 
patterns involve a competition between long- and short- 
range interactions [10] and most theoretical work concen- 
trate on specific examples of such interactions lH-[l5|. 
e.g. spin models with Hamiltonians on the form 



pairwise interactions 



h = kT 



(hj) 



(i) 



^3 



where the spins typically represent some coarse-grained 
feature of the system of interest, for example local charge 
density in a Mott insulator [3] or phases in a Langmuir 
film [16j- However, the qualitative success of such mod- 
els may have little to do with the underlying physics as 
noted by Zaanen in the context of Mott insulators (l7| . 
Indeed, the same general behavior can be observed in 
models with for example only short-ranged, purely re- 
pulsive forces [18]. An explanation for the universality 
of striped morphologies must therefore be independent 
of specific details of the involved forces. The aim of this 
Letter is to present such a general treatment. As ex- 
pected our method shows that stripes appear naturally 
for large classes of models, but the added generality also 
leads to new tools allowing us to design potentials with 
desired properties. 

Here we study a generic Hamiltonian with isotropic 
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where Vij = V(\i — j\) is a matrix representation of the 
potential that only depends on the distance between 
spins i and j, here denoted \i — j\, with (pQ) as a spe- 
cial case. Depending on considerations regarding exper- 
imental fit or theoretical ease, one may take the spins 
in ([2]) to assume continuous values with the restriction 
^isf = TV, corresponding to a spherical model [l9|, or 
take values from some finite set, where s G {±1} and 
s G {0,±1} are the most common choices, equivalent to 
different Ising models [l2[ flBj. 

Consider now the spherical model. Due to the 
quadratic form of the Hamiltonian (|2j) and the constraint, 
s T Vs and s T s = N in matrix notation, the ground state 
is given by an eigenvector of the interaction matrix V cor- 
responding to the lowest (energy) eigenvalue [20]. The 
central observation for our analysis is the existence of 
a common basis of eigenvectors for all radial potentials, 
namely the Fourier basis. To prove this we start by re- 
calling the fact that if two matrices commute, then it is 
possible to find a set of eigenvectors that simultaneously 
diagonalize them [21]. Consider the commutator for two 
interaction matrices V and W: 
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Each term, Viy Wk>j , in the first sum can be represented 
by a triangle, A in Fig. [TJ Assuming that the lattice is 




FIG. 1. Reflection symmetry of triangles on a lattice causes 
the interaction matrices from any two radial potentials to 
commute. 



periodic or infinite, there will for each such triangle exist 
a unique triangle B, constructed as a reflection of A in 
the line equidistant from point i and j (dashed in the 
figure), corresponding to the term WikVkj in the second 
sum. From the reflection symmetry and the pure radial 
dependence of the potential it follows that ViyWk'j — 
WikVkj = 0, which proves that V and W commute. 

It now suffices to find a set of eigenvectors for a partic- 
ular potential. Perhaps the simplest choice is a nearest 
neighbor interaction, Vij = 1 if \i — j\ = 1 and Vij = if 
\i — j\ > 1. If we in addition make an appropriate choice 
of self-interaction, which only shifts the eigenvalues and 
do not affect the eigenvectors, Vu = —2d where d is the 
lattice dimension, Vij becomes a discrete finite difference 
Laplacian on the lattice. It is well known that both the 
discrete and continuous Laplacian have harmonic eigen- 
functions, e.g. f^(x) = Yl^ cos (2nkiXi/L + fa) which 
is an orthogonal eigenbasis in d dimensions when k goes 
over all distances on the reciprocal lattice, L is the linear 
size of the lattice, fa = ±7r/4 and C% is an appropri- 
ate normalization constant. We have thus shown that all 
interaction matrices have a Fourier eigenbasis. An alter- 
native, more direct but for our purposes less illustrative, 
argument for the common Fourier basis is to note that 
the structure of V implies that it is a so called circu- 
lant matrix [22| . for which the result is known in the 
signal processing literature. That the Fourier base effec- 
tively diagonalize the Hamiltonian in the spherical model 
with translationally invariant interactions has also been 
pointed out by Nussinov (23| . 

This result has two important consequences. First, it 
helps us to understand why systems with different in- 
teractions are expected to have similar ground states. 
Second, knowledge of the universal eigenbasis allows 
us to compute the energy spectrum for any particu- 
lar system using a linear transform of the potential. 
From (|2]) it follows that the energy per spin for a har- 
monic eigenfunction with wave vector k are given by 

E $) = wT,x,y V (\^-y\)fk(^)fk(y)' Usin S various 
trigonometric identities and the radial structure of V this 

expression can be reduced to a radial Fourier transform 

d 

E(k) = Y j y{\A)\\^{2i,kin/L) (4) 

r i=l 

where the sum goes over all distances r on the lattice. 
Note that the energy of a configuration can be computed 
through a fast Fourier transform over the lattice. 

The ground state of the spherical model is the eigen- 
vector ft corresponding to the minimum of E(k). The 
simplest ground state patterns in two dimensions are 
checkerboards and stripes with the corresponding wave- 
length, exemplified in Fig. [2^ and c. Further, the sub- 
space of the eigenbasis corresponding to the minimum 
can contain two kinds of degeneracies. First, any change 
of the phases fa leaves the energy invariant. In two 
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FIG. 2. Examples of eigenmodes of the two dimensional 
spherical model, a-c, Any phase shift of a ground state is a 
new ground state, so checkerboards, stripes and everything 
between can be produced by the same model. Linear com- 
binations of /(2,2) with different phase shifts are shown, all 
having the same energy, d, Exchanging the elements of k 
gives a new ground state and linear combinations of them 
give rise to complex morphologies. Shown is |/(3,4) + ^/(4,3)- 



dimensions this means that anything between checker- 
boards and stripes can be produced, as illustrated in 
Fig. EJi-c. Second, the energy is similarly unaffected 
by arbitrary permutations of the elements of /c, reflect- 
ing that the energy only depends on the magnitude of 
the wave vector (seen most clearly in the continuous 
limit (|5])). Linear combinations of vectors with different 
permutations give rise to complex morphologies, exem- 
plified in Fig. EH. 

The eigenmode analysis is exact for the spherical model 
but also has implications for the discrete Ising models. 
It is not directly applicable as in general an eigenvec- 
tor of the interaction matrix cannot be constructed in 
the restricted discrete space of the Ising spins. However, 
continuous eigenvectors are often used to approximate 
solutions to discrete optimization problems, for example 
graph coloring [25j and partitioning networks into mod- 
ules with minimal intra-connectivity 0, [27j- Here we 
use the same strategy to predict ground states for Ising 
spin- 1/2 models with corresponding potentials by map- 
ping the spins in the spherical model to —1 or +1 de- 
pending on their sign: f%(x) = sign(/^(x)). [28] The dis- 
cretization breaks the energy degeneracy and stripes be- 
come energetically favorable compared to checkerboards 
and more complex patterns. To see why we note that in 
each group of degenerate eigenmodes, with wavelength 
|fc|, there exist linear combinations that produce stripes, 
for example cos(fc ■ x). The error introduced by the dis- 
cretization, \\f^(x) — f%{x)\\2 with the standard L 2 norm, 
always increases the energy in the discrete configurations 
when compared to the continuous ground state. Due to 
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FIG. 3. Predicting ground states in different two dimensional Ising models. (Top) Potentials (inset) and their energy spectra 
in | k | -space from the transform ©. a, a purely repulsive potential [18| |, b, two competing interactions |12|] (see equation QJ), 
and c, an RKKY-like interaction 24]. For small k (long wavelengths) the spectra only depend on the magnitude of fc, but for 
large k (short wavelengths) lattice effects breaks the independence on the direction of k and E(\k\) becomes multivalued. To 
illustrate this, we link series with constant k y with lines and in a examples of such states are shown. (Below) Local minima as 
arrived by through Monte Carlo annealing as well as ground states of Ising spin- 1/2 models with corresponding potentials. 



the it-symmetry of the harmonic functions, the differ- 
ence between f^(x) and f^(x) (with appropriate scaling) 
is largest in regions where the continuous function is close 
to 0, i.e. at the interface between + and — regions. From 
this argument it follows that the errortends to be smallest 
for the striped eigenmode since the interface is minimized 
(assuming that the width of the stripes is large compared 
to the lattice spacing). 

There are two exceptions when the ground state does 
not have stripes. Energy spectra with minimum at the 
boundary produce ground states that are either a uni- 
form ferromagnet (the zero frequency mode) or a checker- 
board pattern (the highest frequency mode allowed on 
the lattice) associated with an anti- ferromagnet. These 
two cases can be viewed as degenerate cases of stripes 
with infinite respective infinitesimal width. 

In Fig. [3] some examples of Ising spin- 1/2 models with 
different potentials are shown together with their energy 
spectra in |fc|-space, examples of local minima [H] and 
their ground states. First is a short-ranged, purely repul- 
sive potential related to the model studied in [18]. Sec- 
ond is a nearest neighbor ferromagnet with long-range 
repulsive Coulomb interaction on the form (pQ) from [12| . 
Last is an attenuated Bessel function, Jo(r)/(r + 1), cho- 
sen for its similarity to the RKKY interaction in spin 
glasses 0. We see that, while the potentials are quali- 
tatively very different, the ground states are defined only 
by the minima in the energy spectrum, i.e. by a single 
length scale. Through rescaling, the potentials can be 
adjusted to have identical ground states. This illustrates 
how little observing striped behavior tells us about the 
interactions in a system. The local minima do however 
show a qualitative difference between the potentials in a 



and b and the RKKY-like potential in c, probably related 
to the difference in localization in energy space. 

Equation also has implications for molecular self- 
assembly. The Fourier basis in the transformation is or- 
thogonal and can be inverted to find the potential corre- 
sponding to a given energy spectrum. This allows us to 
design, from an observed striped state, families of poten- 
tials that generate similar patterns at low-temperature 
by identifying the dominant wavelength and invert an 
energy spectrum with a minimum at this wavelength. 
A demonstration of the procedure is shown in Fig. 2J 
the Fourier power spectrum of a pixelised image of a 
metastable state in an experimental system [6] was calcu- 
lated; an energy spectrum was constructed with gaussian 
minimum at the same wavelength as the experimental 
system; and finally the corresponding potential was found 
using the inverse transform of (j4]). The constructed sys- 
tem has striped metastable states similar to those found 
in the experimental system. We conclude that it is rela- 
tively easy to construct families of potentials with desired 
metastable striped morphologies. 

In the continuous limit the transformation ((4]) becomes 
a Hankel transform, in two dimensions defined as 

PCX) 

E(k)=27r rdrV(r)J (27r\k\r) (5) 
Jo 

where Jo is a Bessel function of the first kind. For the 
general expression in higher dimensions, see [30j ]. Note 
that in the continuum limit the energy only depends 
on the magnitude of the wavevector since the effects of 
the principal lattice directions disappear. As noted in 
Fig. [31 this independence holds true on the lattice as 
well for small wavevectors and sufficiently long-range in- 






FIG. 4. Designing interactions to imitate observed striped patterns, a, Stripes in a ferrofluid confined between two glass plates 
in a magnetic field [from 0, @], reprinted with permission from A A AS and Elsevier] b, the (negative) radial power spectrum 
(blue) of previous picture together with a gaussian (red), (inset) a pair potential V(r) constructed as the (inverse) transform J4| 
of said gaussian and c, metastable state of an Ising spin-1/2 model with potential V(r). Note that the chosen energy spectrum is 
not unique. Many potentials having a spectrum minimized at the same wavelength will show similar low-temperature behavior. 



ter act ions. Equation (|5j) allows us to use the analytical 
properties of the Hankel transform to qualitatively un- 
derstand for example why the Bessel function of Fig. [3t 
has such a sharp spectrum: the Hankel transform of a 
Bessel function is a Dirac delta function. 

In summary we have shown that the energy spectrum 
of spherical spin systems with isotropic interactions can 
be derived directly from the Fourier transform of the po- 
tential. Due to a degeneracy in the energy eigenstates 
the spherical model has ground states with various pat- 
terns such as stripes, checkerboards, and more compli- 
cated morphologies. In discrete spin models the degener- 
acy is broken leading to striped ground states being ener- 
getically favored. We suggest that this can offer a generic 
explanation to why striped patterns are so frequently ob- 
served in various experimental and natural systems. 

The authors would like to thank Olle Haggstrom for 
pointing out how purely repulsive potentials can give rise 
to striped ground states. 
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